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ABSTRACT 

The nearby dwarf spheroidal galaxy Draco with its high mass to light ratio 
is one of the most auspicious targets for indirect dark matter (DM) searches. 
Annihilation of hypothetical DM particles can result in high-energy 7-rays, e.g. 
from neutralino annihilation in the supersymmetric framework. 
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With the MAGIC telescope a search for a possible DM signal originating from 
Draco was performed during 2007. The analysis of the data results in a flux upper 
limit (2cr) of 1.1 x 10~^^ photons cm~^ sec~^ for photon energies above 140 GeV, 
assuming a point like source. Furthermore, a comparison with predictions from 
supersymmetric models is given. While our results do not constrain the mSUGRA 
phase parameter space, a very high flux enhancement can be ruled out. 

Subject headings: dark matter — gamma rays: observations — galaxies: dwarf 
— galaxies: individual (Draco) 



Introduction 



Astronomical observations provide strong evidence for the existence of a new type of 
non-luminous, non-baryonic matter, contrib uting to the total eii ergy density of the universe 
about six times more than baryonic matter (jSpergel et al.l 120071 ). This so-called Dark Mat- 
ter (DM) makes its presence known through gravitational effects, and could be made of 
so-far undetected relic particles from the Big Bang. Weakly Interacting Massive Particles 
(WIMP) are candidates for DM, with the lightest supersymmetric particle (neutralino) be- 
ing one of the most favored candidates in the list of possible WIMPs. Stable neutralinos 
are predi cted in many supersymmetr i c (SU SY) extensions of the Standard Model of Particle 
Physics ( iJungman fc Kamionkowskil Il995l ). Since the neutralinos are Majorana particles, 
pairs of neutralinos can annihilate and produce Standard Model particles. Direct annihila- 
tions into 77 or Z7 produce a sharp line spectrum with a photon energy depending on the 
neutralino mass. Unfortunately, these processes are loop-suppressed and therefore very rare. 
Neutralinos can also annihilate to pairs of r or quarks, leading in subsequent processes to 
7r°-decays, resulting in a continuous photon spectrum. 

Draco is a dwarf spheroidal galaxy accompanying the Milky Way at a galactocentric distance 
of about 82kpc. It is ch aracterized by a h i gh mass to light ratio M/L > 200, implying a 
high DM concentration (IMayer et al.l 120071 : iBergstrom fc Hooperl 120061) so corn plying with 
the trend generally deduced for low-luminosity galaxies (e.g. iPersic et al.l (119961 )). 



2. Expected 7-Ray Flux From Neutralino Self- Annihilation 

The expected 7-ray flux depends on details of the supersymmetric (SUSY) model as 
well as on the density distribution of the DM in the observed source. In general, the DM 
is assumed to be distributed in an extended halo around spheroidal galaxies. The radial 
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profile of the DM distribution in the halo is modeled by a power law, pBuir) = Cr~'^, where 
the parameter e > describes the shape of the DM distribution in the crucial innermost 
region, e = results in a so-called core model with a central flat region, whereas profiles 
with 0.7 < e < 1.2 der iote the so-called c usp p rofiles. In addition, we chose an exponential 
cut-off as proposed by iKazantzidis et al.l (120041 ): 



Pdm 



Cr' 



exp 



r 



with the values for r?,, C and e given in table 1 for a cusp and a core profile for Draco. 
With the present angular resolution of the MAGIC telescope (0.1°), the two models are 
indistinguishable due to the limited angular resolution of the telescope which smears the 
determination of the profile (see figure [1]). From this figure we can see that the two pro- 
files are discriminated at an angular distance of 0.4°, where the intrinsic flux is already 
decreased by a factor 20. Depending on SUSY model parameters, the annihilation cross 
section, the average number of photons produced per annihilation and the shape of the 7 
spectrum can drastically change. Also a change in the shape of the DM density distribu- 
tion along the line of sight can significantly change the 7-ray flux. Formula (1) describes 
the expected 7-ray flux above an energy Eq from neutralino self-annihilation within Draco. 



where : 

/busy 

p(r) 

(av) 
los 



^/SUSY < ^(^) >An 
N~f{E>Eo){av) 

DM density profile derived for Draco 
photon yield per annihilation with E > Eq 
neutralino mass 

thermally averaged annihilation cross-section 
Point Spread Function (PSF) of the telescope 
Pointing angle. (\E' = for the center of Draco) 
solid angle of the telescope's resolution 
the line of sight 



The factor < Ji"^) >n is shown for the cusp and the core model in figure [H Even 
though this factor converges for both models for small pointing angles to the same value 
< J(0) >n, there is an uncertainty in the distribution of the DM b y the existence of a 
hypothetical central blac k hole or a clumpy distribution of the DM (jStrigari et al.l 120071 : 
Colafrancesco et al.l 120071 ). which could lead to a significant flux enhancement. 
Due to the high predictive power of the mSUGRA framework, where the SUSY breaking 
effect s are transmitted from the high energy sca le to th e electroweak scale by the gravi- 
ton, Jchamseddine et aLlllQsi llnoue et al.l[l982al jbl. Il984h we simulated several million mod- 
els using the following parameters: mo < 6 TeV, mi/2 < 4 TeV, -4 TeV < Aq < 4 TeV, 
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tan/? < 50 and /i > (IStark et al.ll2005l : iGondolo et al.ll2000l ). Figure [2] summarizes the re- 
sulting thermally averaged neutralino annihilation cross sections for all models not violating 
any observational constraints as well as resulting in a total DM relic density I^dm^^ in agree- 
ment wit h the 2(7 upper limit (u.l.) of 0.113 as derived from combined data from SPSS and 
WMAP dXeemark et aPbood ). The yellow points correspond to models with ttiq < 2 TeV 
(as favored by particle physics), and the blue points represent mo > 2 TeV. Models resulting 
in a relic density below the lower WMAP-limit of 0.097 are included, since neutralinos could 
contribute only a fraction to the total DM in the universe. For these models (shown as dark 

blue points and dark yellow points in figure [2]), a scale factor of n = is applied 

to adjust for the DM relic density. 



3. Observation of Draco and Analysis 



Among all Imagin g Air Cherenkoy Tele s copes in operation, M AGIC is the largest single- 
dish facility (see e.g. (IBaixeras et al.ll2004l : ICortina et al.ll2005l ) for a detailed description) 
and has the lowest energy threshold. MAGIC is located on the Canary Island La Palma 
(28.8°N, 17.8°W, 2200 m a.s.L). The field of view (FOV) of the 576-pixel photomultiplier 
camera is 3.5°. The angular resolution is ~ 0.1° and the energy resolution above 150 GeV is 
about 25%. MAGIC has a trigger threshold of ~60 GeV for small zenith angles (ZA), which 
increases for larger ZAs. 



Data were taken in the false-source tracking (wobble) mode (IFomin et al.l 1 19941 ) with two 
pointing directions at 24' distance and opposite sides of the source direction in May 2007 for 
a total observation time of 7.8 hours. Even though the source is expected to be extended, 
the wobble mode is justified, as at a distance of 24' of the center of Draco the expected flux 
from this direction is less than 5% of the flux coming from the center of Draco for both the 
cusp and the core model (see flgure[T]). The ZA ranges between 29° and 42°. 



Firstly, the calibration of the data ( Gaug|[2"o05l ) was performed. The arrival times of the 
photons in core pixels (> 6 photoelectrons (phe)) are required to be within a time window 
of 4.5 ns and for boundary pixels (> 3 phe) within a time window of 1.5 ns of a neigh boring 
core pixel. The next step includes the Hillas parameterization of the shower images (IHillas 
19851 ). Two additional parameters, namely the time gradi ent along the main shower axis 
and the time spread of the shower pixels, were computed (ITescaro et al.l 120071) . Hadronic 



background suppre ssion was achieved using the Random Forest (RF) method (jBreimanll2001 



Albert et al.l 120071 ). where for each event the so-called Hadronness is computed, based on 
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the Hillas and the time parameters. 7/hadron separatio n is reahzed by a cut in Hadron- 



NESS, derived from a 7 Monte Carlo (MC) test sample (IHeck et al.lll998l : iMajumdar et al. 



2OO5I ) . requiring a 7-cut efficiency of 70%. Moreover, the RF method was also used for the 



energy estimation. 



4. Results 



We searched for a steady 7-ray emission from the direction of the dwarf spheroidal 
galaxy Draco. The analysis energy threshold defined as the peak of the energy distribution 
of MC generated 7 events after cuts is 140 GeV. Images of air showers initiated by 7-rays 
coming from the center of Draco are characterized by a small a image parameter, which is 
the angle between the main axis of the shower image and the connecting line between the 
center of gravity of the shower and the source position in the camera. The distribution of a 
is shown in figure [3] for all events after cuts. No significant exce ss was found . The 2a u.l. on 
the number of excess events was calculated using the method of iRolke et al.l (120051 ) . applying 
a systematic error of 30%. The number of excess events were converted into an integral flux 
U.L, depending on the assumed underlying spectrum. For a power law with spectral index 
—1.5, typical for a DM annihilation spectrum, and assuming a point like source, the 2a u.l. 
is: 

^2a{E > 140 GeV) = 1.1 X 10"^^ photons cm"^ s~\ 



For different mSUGRA model parameters using the benchmark points defined in ( iBattaglia et al 



20041 ) and for other models, we computed the 7-ray spectra expected from neutralino anni- 
hilations. Assuming these underlying spectra, the u.l. on the integrated flux above 140 GeV 
is computed. Using formula (1) and assuming a DM distribution following the profiles ac- 
cording to table 1, the flux u.l. is displayed in the units of a thermally averaged cross section 
in table 2 and in figure [2j 



'S>^(E>Eo) 



< av > 



As can be seen, the measured flux u.l. is several orders of magnitude larger than pre- 
dicted for the smooth DM distribution. But a high clumpy structure of the DM distribution 
or a central black hole c ould provide a significant flux enhancement ( jStrigari et al.l 120071 : 
Colafrancesco et al.l 120071 ). which would decrease ^2^. The analysis presented here can set 
a limit on the flux enhancement depending on the mSUGRA input parameters. For the 
benchmark models the values for k < av > and are displayed in table 2 and in figure [2j 
For these models, the u.l. on the flux enhancement is around 0(10^ — 10^). 
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Conclusions 



We present the first search for 7-rays from the direction of Draco using an Imaging 
Air Cherenkov Telescope (lACT). No signal was detected. The 2o" u.l. on a steady 7-ray 
emission above 140 GeV originating from Draco does not exceed 1.1 x 10^^^ photons cm~^ s~^ 
if the underlying spectrum follows a power law with spectral index - 1.5. 
For the mSUGRA benchmark models defined in (IBattaglia et al.l l2004r) and as suming a 
smooth DM density distribution for Draco as given in (jSanchez-Conde et al.ll2007l ). our flux 
upper hmits are 0(10^ - 10^) above the predicted values. It is therefore not possible to 
constrain the mSUGRA phase space by these results, but a very high flux enhancement can 
be excluded. 

Even though an indirect DM detection by measuring 7-rays from neutralino annihilation 
within the halo of Draco seems for present lACTs out of reach, future satellite telescopes 
like GLAST with lower energy thresholds might be sensitive enough to reach the mSUGRA 
parameter space. 
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profile 


e 




C 


r6(kpc) 


cusp 


1 


3.1 > 


< lO'^ Mq kpc"2 


1.189 


core 





3.6 > 


< 10^ Mq kpc"^ 


0.238 



Table 1: Parameters considered for cusp and core DM density profiles (jSanchez-Conde et al, 
2007h . 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 
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model: 


A' 


B' 


C' 


D' 


E' 


F' 


G' 


H' 


mo [GcV] 


107 


57 


80 


101 


1532 


3440 


113 


244 


Ao [GcV] 


























nii/2 [GcV] 


600 


250 


400 


525 


300 


1000 


375 


935 


tan f3 


5 


10 


10 


10 


10 


10 


20 


20 


[GcV] 


243 


95 


158 


212 


121 


428 


148 


389 


K < (7V > 

[cm^s-1] 


5.55 • 10"^^ 


6.83 • 10"^** 


1.42 ■ 10"^** 


6.05 ■ 10"^^ 


3.74 ■ 10"^" 


3.65 • 10"^° 


7.18 • 10"^** 


2.87 • lO"^** 




















1.63 • 10"^^ 


1.79 • 10"^^ 


1.00 ■ 10"^^ 


1.15 ■ 10"^^ 


4.82 ■ 10"^^ 


1.42 • 10"^^ 


1.11 • 10"^^ 


3.92 • 10"^^ 


[cm^s-i] 


















ul on fiux 
enhancement 


2.9 • 10^ 


2.6 • lO" 


7.0 ■ lO'^ 


1.9 ■ 10^ 


1.3 ■ lO" 


3.9 • 10^ 


1.5 • 10^ 


1.4 ■ lO** 


model: 


r 


J' 


K' 


L' 


a 


/3 


7 


5 


mo [GcV] 


181 


299 


1001 


303 


5980 


180 


1140 


4540 


Ao [GcV] 














-300 


-2800 


-1800 


300 


mi/2 [GcV] 


350 


750 


1300 


450 


680 


720 


1120 


300 


tan P 


35 


35 


46 


47 


50 


5 


50 


35 


[GcV] 


138 


309 


554 


181 


277 


301 


479 


109 


At < (7V > 

[cm^s-i] 


3.17 • 10"^'' 


1.67 • 10"^** 


2.22 ■ 10"^'' 


3.85 ■ 10"^^ 


2.27 ■ lO"^** 


9.75 • 10"^** 


9.43 • 10"^^ 


1.79 • 10"^* 


















[em^s-i] 


1.08 • 10"^"^ 


9.02 • 10"^^ 


1.44 ■ 10"^^ 


1.91 ■ 10"^^ 


7.22 ■ 10"^^ 


4.78 • 10"^^ 


2.97 • 10"^^ 


4.39 • 10"^° 


















ul on flux 
enhancement 


3.4 • 10^ 


5.4 • 10^ 


6.5 ■ 10^ 


4.9 ■ lO'^ 


3.2 ■ lO*' 


4.9 • lO'"" 


3.1 • lO'"" 


2.4 ■ 10*^ 



Table 2: Thermally averaged neutralino annihilation cross section < av >, the u.l. on the 
flux F2a, displayed in units of < av >, and the 2a u.l. on the flux enhan cement. Models A' 
- L' correspond to the benchmark models given by (iBattaglia et al.ll2004j ). Models a - 6 are 
typical models chosen by the authors with Aq ^ 0. 
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Fig. 1. — The factor < Ji^^) >An for the cusp (red) and the core (blue) profile. \1' = 
corresponds to the center of Draco. At an angular distance of 0.4° of the center of Draco, < 
</(^) >An is reduced by a factor of around 20 for both models. 
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Fig. 2. — Thermally averaged neutralino annihilation cross section as a function of the 
neutralino mass for mSUGRA models after renormalization to the relic density, as de- 
scri bed in the text. The red dots marked with roman letters indicate benchmark models 
by (IBattaglia et al.ll2004l ). The ones marked with greek letters are models chosen by the 
authors. The red boxes indicate the flux upper hmit, displayed in units of < crw >, assuming 
a smooth DM halo as given in table 1. See text for details. 
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Fig. 3. — Distribution of the a parameter for 7-ray candidates coming from the center of 
Draco (black marker) and background (grey markers) for data taken between 05/09/2007 
- 05/20/2007. The energy threshold is 140 GeV. For the signal region [a < 12°), are the 
number of ON events: 10883 and the number of OFF events: 10996. With the method of 



( iRolke et al.ll2005l ). the 2a upper limit on excess events is 231. 



